RNA three-dimensional structure drives the sequence organization of potato spindle tuber viroid quasispecies.

RNA viruses and viroids exist and evolve as quasispecies due to error-prone replication. Quasispecies consist of a few dominant master sequences alongside numerous variants that contribute to genetic diversity. Upon environmental changes, certain variants within quasispecies have the potential to become the dominant sequences, leading to the emergence of novel infectious strains. However, the emergence of new infectious variants remains unpredictable. Using mutant pools prepared by saturation mutagenesis of selected stem and loop regions, our study of potato spindle tuber viroid (PSTVd) demonstrates that mutants forming local three-dimensional (3D) structures similar to the wild type (WT) are more likely to accumulate in PSTVd quasispecies. The selection mechanisms underlying this biased accumulation are likely associated with cell-to-cell movement and long-distance trafficking. Moreover, certain trafficking-defective PSTVd mutants can be spread by functional sister genomes in the quasispecies. Our study reveals that the RNA 3D structure of stems and loops constrains the evolution of viroid quasispecies. Mutants with a structure similar to WT have a higher likelihood of being maintained within the quasispecies and can potentially give rise to novel infectious variants. These findings emphasize the potential of targeting RNA 3D structure as a more robust approach to defend against viroid infections.

Unveiling hidden structural patterns in the SARS-CoV-2 genome: Computational insights and comparative analysis.

SARS-CoV-2, the causative agent of COVID-19, is known to exhibit secondary structures in its 5' and 3' untranslated regions, along with the frameshifting stimulatory element situated between ORF1a and 1b. To identify additional regions containing conserved structures, we utilized a multiple sequence alignment with related coronaviruses as a starting point. We applied a computational pipeline developed for identifying non-coding RNA elements. Our pipeline employed three different RNA structural prediction approaches. We identified forty genomic regions likely to harbor structures, with ten of them showing three-way consensus substructure predictions among our predictive utilities. We conducted intracomparisons of the predictive utilities within the pipeline and intercomparisons with four previously published SARS-CoV-2 structural datasets. While there was limited agreement on the precise structure, different approaches seemed to converge on regions likely to contain structures in the viral genome. By comparing and combining various computational approaches, we can predict regions most likely to form structures, as well as a probable structure or ensemble of structures. These predictions can be used to guide surveillance, prophylactic measures, or therapeutic efforts. Data and scripts employed in this study may be found at https://doi.org/10.5281/zenodo.8298680.

Conserved structures and dynamics in 5'-proximal regions of Betacoronavirus RNA genomes.

Betacoronaviruses are a genus within the Coronaviridae family of RNA viruses. They are capable of infecting vertebrates and causing epidemics as well as global pandemics in humans. Mitigating the threat posed by Betacoronaviruses requires an understanding of their molecular diversity. The development of novel antivirals hinges on understanding the key regulatory elements within the viral RNA genomes, in particular the 5'-proximal region, which is pivotal for viral protein synthesis. Using a combination of cryo-electron microscopy, atomic force microscopy, chemical probing, and computational modeling, we determined the structures of 5'-proximal regions in RNA genomes of Betacoronaviruses from four subgenera: OC43-CoV, SARS-CoV-2, MERS-CoV, and Rousettus bat-CoV. We obtained cryo-electron microscopy maps and determined atomic-resolution models for the stem-loop-5 (SL5) region at the translation start site and found that despite low sequence similarity and variable length of the helical elements it exhibits a remarkable structural conservation. Atomic force microscopy imaging revealed a common domain organization and a dynamic arrangement of structural elements connected with flexible linkers across all four Betacoronavirus subgenera. Together, these results reveal common features of a critical regulatory region shared between different Betacoronavirus RNA genomes, which may allow targeting of these RNAs by broad-spectrum antiviral therapeutics.

The application of single-molecule optical tweezers to study disease-related structural dynamics in RNA.

RNA, a dynamic and flexible molecule with intricate three-dimensional structures, has myriad functions in disease development. Traditional methods, such as X-ray crystallography and nuclear magnetic resonance, face limitations in capturing real-time, single-molecule dynamics crucial for understanding RNA function. This review explores the transformative potential of single-molecule force spectroscopy using optical tweezers, showcasing its capability to directly probe time-dependent structural rearrangements of individual RNA molecules. Optical tweezers offer versatility in exploring diverse conditions, with the potential to provide insights into how environmental changes, ligands and RNA-binding proteins impact RNA behaviour. By enabling real-time observations of large-scale structural dynamics, optical tweezers emerge as an invaluable tool for advancing our comprehension of RNA structure and function. Here, we showcase their application in elucidating the dynamics of RNA elements in virology, such as the pseudoknot governing ribosomal frameshifting in SARS-CoV-2.

Bipartite viral RNA genome heterodimerization influences genome packaging and virion thermostability.

Multi-segmented viruses often multimerize their genomic segments to ensure efficient and stoichiometric packaging of the correct genetic cargo. In the bipartite Nodaviridae family, genome heterodimerization is also observed and conserved among different species. However, the nucleotide composition and biological function for this heterodimer remain unclear. Using Flock House virus as a model system, we developed a next-generation sequencing approach ("XL-ClickSeq") to probe heterodimer site sequences. We identified an intermolecular base-pairing site which contributed to heterodimerization in both wild-type and defective virus particles. Mutagenic disruption of this heterodimer site exhibited significant deficiencies in genome packaging and encapsidation specificity to viral genomic RNAs. Furthermore, the disruption of this intermolecular interaction directly impacts the thermostability of the mature virions. These results demonstrate that the intermolecular RNA-RNA interactions within the encapsidated genome of an RNA virus have an important role on virus particle integrity and thus may impact its transmission to a new host.IMPORTANCEFlock House virus is a member of Nodaviridae family of viruses, which provides a well-studied model virus for non-enveloped RNA virus assembly, cell entry, and replication. The Flock House virus genome consists of two separate RNA molecules, which can form a heterodimer upon heating of virus particles. Although similar RNA dimerization is utilized by other viruses (such as retroviruses) as a packaging mechanism and is conserved among Nodaviruses, the role of heterodimerization in the Nodavirus replication cycle is unclear. In this research, we identified the RNA sequences contributing to Flock House virus genome heterodimerization and discovered that such RNA-RNA interaction plays an essential role in virus packaging efficiency and particle integrity. This provides significant insight into how the interaction of packaged viral RNA may have a broader impact on the structural and functional properties of virus particles.

The molecular language of RNA 5' ends: guardians of RNA identity and immunity.

RNA caps are deposited at the 5' end of RNA polymerase II transcripts. This modification regulates several steps of gene expression, in addition to marking transcripts as self to enable the innate immune system to distinguish them from uncapped foreign RNAs, including those derived from viruses. Specialized immune sensors, such as RIG-I and IFITs, trigger antiviral responses upon recognition of uncapped cytoplasmic transcripts. Interestingly, uncapped transcripts can also be produced by mammalian hosts. For instance, 5'-triphosphate RNAs are generated by RNA polymerase III transcription, including tRNAs, Alu RNAs, or vault RNAs. These RNAs have emerged as key players of innate immunity, as they can be recognized by the antiviral sensors. Mechanisms that regulate the presence of 5'-triphosphates, such as 5'-end dephosphorylation or RNA editing, prevent immune recognition of endogenous RNAs and excessive inflammation. Here, we provide a comprehensive overview of the complexity of RNA cap structures and 5'-triphosphate RNAs, highlighting their roles in transcript identity, immune surveillance, and disease.

Diversity and modularity of tyrosine-accepting tRNA-like structures.

Certain positive-sense single-stranded RNA viruses contain elements at their 3' termini that structurally mimic tRNAs. These tRNA-like structures (TLSs) are classified based on which amino acid is covalently added to the 3' end by host aminoacyl-tRNA synthetase. Recently, a cryoEM reconstruction of a representative tyrosine-accepting tRNA-like structure (TLSTyr) from brome mosaic virus (BMV) revealed a unique mode of recognition of the viral anticodon-mimicking domain by tyrosyl-tRNA synthetase. Some viruses in the hordeivirus genus of Virgaviridae are also selectively aminoacylated with tyrosine, yet these TLS RNAs have a different architecture in the 5' domain that comprises the atypical anticodon loop mimic. Herein, we present bioinformatic and biochemical data supporting a distinct secondary structure for the 5' domain of the hordeivirus TLSTyr compared to those in Bromoviridae Despite forming a different secondary structure, the 5' domain is necessary to achieve robust in vitro aminoacylation. Furthermore, a chimeric RNA containing the 5' domain from the BMV TLSTyr and the 3' domain from a hordeivirus TLSTyr are aminoacylated, illustrating modularity in these structured RNA elements. We propose that the structurally distinct 5' domain of the hordeivirus TLSTyrs performs the same role in mimicking the anticodon loop as its counterpart in the BMV TLSTyr Finally, these structurally and phylogenetically divergent types of TLSTyr provide insight into the evolutionary connections between all classes of viral tRNA-like structures.

Structural flexibility of Toscana virus nucleoprotein in the presence of a single-chain camelid antibody.

Phenuiviridae nucleoprotein is the main structural and functional component of the viral cycle, protecting the viral RNA and mediating the essential replication/transcription processes. The nucleoprotein (N) binds the RNA using its globular core and polymerizes through the N-terminus, which is presented as a highly flexible arm, as demonstrated in this article. The nucleoprotein exists in an `open' or a `closed' conformation. In the case of the closed conformation the flexible N-terminal arm folds over the RNA-binding cleft, preventing RNA adsorption. In the open conformation the arm is extended in such a way that both RNA adsorption and N polymerization are possible. In this article, single-crystal X-ray diffraction and small-angle X-ray scattering were used to study the N protein of Toscana virus complexed with a single-chain camelid antibody (VHH) and it is shown that in the presence of the antibody the nucleoprotein is unable to achieve a functional assembly to form a ribonucleoprotein complex.

Conformational transitions of the HIV-1 Gag polyprotein upon multimerization and gRNA binding.

During the HIV-1 assembly process, the Gag polyprotein multimerizes at the producer cell plasma membrane, resulting in the formation of spherical immature virus particles. Gag-genomic RNA (gRNA) interactions play a crucial role in the multimerization process, which is yet to be fully understood. We performed large-scale all-atom molecular dynamics simulations of membrane-bound full-length Gag dimer, hexamer, and 18-mer. The inter-domain dynamic correlation of Gag, quantified by the heterogeneous elastic network model applied to the simulated trajectories, is observed to be altered by implicit gRNA binding, as well as the multimerization state of the Gag. The lateral dynamics of our simulated membrane-bound Gag proteins, with and without gRNA binding, agree with prior experimental data and help to validate our simulation models and methods. The gRNA binding is observed to affect mainly the SP1 domain of the 18-mer and the matrix-capsid linker domain of the hexamer. In the absence of gRNA binding, the independent dynamical motion of the nucleocapsid domain results in a collapsed state of the dimeric Gag. Unlike stable SP1 helices in the six-helix bundle, without IP6 binding, the SP1 domain undergoes a spontaneous helix-to-coil transition in the dimeric Gag. Together, our findings reveal conformational switches of Gag at different stages of the multimerization process and predict that the gRNA binding reinforces an efficient binding surface of Gag for multimerization, and also regulates the dynamic organization of the local membrane region itself.

The 3' terminal region of Zika virus RNA contains a conserved G-quadruplex and is unfolded by human DDX17.

Zika virus (ZIKV) infection remains a worldwide concern, and currently no effective treatments or vaccines are available. Novel therapeutics are an avenue of interest that could probe viral RNA-human protein communication to stop viral replication. One specific RNA structure, G-quadruplexes (G4s), possess various roles in viruses and all domains of life, including transcription and translation regulation and genome stability, and serves as nucleation points for RNA liquid-liquid phase separation. Previous G4 studies on ZIKV using a quadruplex forming G-rich sequences Mapper located a potential G-quadruplex sequence in the 3' terminal region (TR) and was validated structurally using a 25-mer oligo. It is currently unknown if this structure is conserved and maintained in a large ZIKV RNA transcript and its specific roles in viral replication. Using bioinformatic analysis and biochemical assays, we demonstrate that the ZIKV 3' TR G4 is conserved across all ZIKV isolates and maintains its structure in a 3' TR full-length transcript. We further established the G4 formation using pyridostatin and the BG4 G4-recognizing antibody binding assays. Our study also demonstrates that the human DEAD-box helicases, DDX3X132-607 and DDX17135-555, bind to the 3' TR and that DDX17135-555 unfolds the G4 present in the 3' TR. These findings provide a path forward in potential therapeutic targeting of DDX3X or DDX17's binding to the 3' TR G4 region for novel treatments against ZIKV.

Exploring the conformational dynamics of the SARS-CoV-2 SL4 hairpin by combining optical tweezers and base analogues.

The parasitic nature of the SARS-CoV-2 virus demands selective packaging of its RNA genome (gRNA) from the abundance of other nucleic acids present in infected cells. Despite increasing evidence that stem-loop 4 (SL4) of the gRNA 5' UTR is involved in the initiation of this process by binding the nucleocapsid (N) protein, little is known about its conformational dynamics. Here, we unravel the stability, dynamics and (un)folding pathways of SL4 using optical tweezers and a base analogue, tCO, that provides a local and subtle increase in base stacking without perturbing hydrogen bonding. We find that SL4 (un)folds mainly in a single step or through an intermediate, encompassing nucleotides from the central U bulge to the hairpin loop. Due to an upper-stem CU mismatch, SL4 is prone to misfold, the extent of which can be tuned by incorporating tCO at different positions. Our study contributes to a better understanding of SARS-CoV-2 packaging and the design of drugs targeting SL4. We also highlight the generalizability of using base analogues in optical tweezers experiments for probing intramolecular states and conformational transitions of various nucleic acids at the level of single molecules and with base-pair resolution.

Structural and computational studies of HIV-1 RNA.

Viruses remain a global threat to animals, plants, and humans. The type 1 human immunodeficiency virus (HIV-1) is a member of the retrovirus family and carries an RNA genome, which is reverse transcribed into viral DNA and further integrated into the host-cell DNA for viral replication and proliferation. The RNA structures from the HIV-1 genome provide valuable insights into the mechanisms underlying the viral replication cycle. Moreover, these structures serve as models for designing novel therapeutic approaches. Here, we review structural data on RNA from the HIV-1 genome as well as computational studies based on these structural data. The review is organized according to the type of structured RNA element which contributes to different steps in the viral replication cycle. This is followed by an overview of the HIV-1 transactivation response element (TAR) RNA as a model system for understanding dynamics and interactions in the viral RNA systems. The review concludes with a description of computational studies, highlighting the impact of biomolecular simulations in elucidating the mechanistic details of various steps in the HIV-1's replication cycle.

A novel tamanavirus (Flaviviridae) of the European common frog (Rana temporaria) from the UK.

Flavivirids are small, enveloped, positive-sense RNA viruses from the family Flaviviridae with genomes of ~9-13 kb. Metatranscriptomic analyses of metazoan organisms have revealed a diversity of flavivirus-like or flavivirid viral sequences in fish and marine invertebrate groups. However, no flavivirus-like virus has been identified in amphibians. To remedy this, we investigated the virome of the European common frog (Rana temporaria) in the UK, utilizing high-throughput sequencing at six catch locations. De novo assembly revealed a coding-complete virus contig of a novel flavivirid ~11.2 kb in length. The virus encodes a single ORF of 3456 aa and 5' and 3' untranslated regions (UTRs) of 227 and 666 nt, respectively. We named this virus Rana tamanavirus (RaTV), as BLASTp analysis of the polyprotein showed the closest relationships to Tamana bat virus (TABV) and Cyclopterus lumpus virus from Pteronotus parnellii and Cyclopterus lumpus, respectively. Phylogenetic analysis of the RaTV polyprotein compared to Flavivirus and Flavivirus-like members indicated that RaTV was sufficiently divergent and basal to the vertebrate Tamanavirus clade. In addition to the Mitcham strain, partial but divergent RaTV, sharing 95.64-97.39 % pairwise nucleotide identity, were also obtained from the Poole and Deal samples, indicating that RaTV is widespread in UK frog samples. Bioinformatic analyses of predicted secondary structures in the 3'UTR of RaTV showed the presence of an exoribonuclease-resistant RNA (xrRNA) structure standard in flaviviruses and TABV. To examine this biochemically, we conducted an in vitro Xrn1 digestion assay showing that RaTV probably forms a functional Xrn1-resistant xrRNA.

An RNA excited conformational state at atomic resolution.

Sparse and short-lived excited RNA conformational states are essential players in cell physiology, disease, and therapeutic development, yet determining their 3D structures remains challenging. Combining mutagenesis, NMR spectroscopy, and computational modeling, we determined the 3D structural ensemble formed by a short-lived (lifetime ~2.1 ms) lowly-populated (~0.4%) conformational state in HIV-1 TAR RNA. Through a strand register shift, the excited conformational state completely remodels the 3D structure of the ground state (RMSD from the ground state = 7.2 ± 0.9 Å), forming a surprisingly more ordered conformational ensemble rich in non-canonical mismatches. The structure impedes the formation of the motifs recognized by Tat and the super elongation complex, explaining why this alternative TAR conformation cannot activate HIV-1 transcription. The ability to determine the 3D structures of fleeting RNA states using the presented methodology holds great promise for our understanding of RNA biology, disease mechanisms, and the development of RNA-targeting therapeutics.

Identification of a putative Gag binding site critical for feline immunodeficiency virus genomic RNA packaging.

The retroviral Gag precursor plays a central role in the selection and packaging of viral genomic RNA (gRNA) by binding to virus-specific packaging signal(s) (psi or ψ). Previously, we mapped the feline immunodeficiency virus (FIV) ψ to two discontinuous regions within the 5' end of the gRNA that assumes a higher order structure harboring several structural motifs. To better define the region and structural elements important for gRNA packaging, we methodically investigated these FIV ψ sequences using genetic, biochemical, and structure-function relationship approaches. Our mutational analysis revealed that the unpaired U85CUG88 stretch within FIV ψ is crucial for gRNA encapsidation into nascent virions. High-throughput selective 2' hydroxyl acylation analyzed by primer extension (hSHAPE) performed on wild type (WT) and mutant FIV ψ sequences, with substitutions in the U85CUG88 stretch, revealed that these mutations had limited structural impact and maintained nucleotides 80-92 unpaired, as in the WT structure. Since these mutations dramatically affected packaging, our data suggest that the single-stranded U85CUG88 sequence is important during FIV RNA packaging. Filter-binding assays performed using purified FIV Pr50Gag on WT and mutant U85CUG88 ψ RNAs led to reduced levels of Pr50Gag binding to mutant U85CUG88 ψ RNAs, indicating that the U85CUG88 stretch is crucial for ψ RNA-Pr50Gag interactions. Delineating sequences important for FIV gRNA encapsidation should enhance our understanding of both gRNA packaging and virion assembly, making them potential targets for novel retroviral therapeutic interventions, as well as the development of FIV-based vectors for human gene therapy.

Identification, Sequencing, and Molecular Analysis of RNA2 of Artichoke Italian Latent Virus Isolates from Known Hosts and a New Host Plant Species.

Despite its first description in 1977 and numerous reports of its presence in various plant species in many countries, the molecular information available in GenBank for artichoke Italian latent virus (AILV) is still limited to a single complete genome sequence (RNA1 and 2) of a grapevine isolate (AILV-V) and a partial portion of the RNA2 sequence from an isolate of unknown origin and host. Here, we report the results of molecular analyses conducted on the RNA2 of some AILV isolates, sequenced for the first time in this study, together with the first-time identification of AILV in a new host plant species, namely chard (Beta vulgaris subsp. vulgaris), associated with vein clearing and mottling symptoms on leaves. The different AILV isolates sequenced were from artichoke (AILV-C), gladiolus (AILV-G), Sonchus (AILV-S), and chard (AILV-B). At the molecular level, the sequencing results of the RNA2 segments showed that AILV-C, AILV-G, AILV-S, and AILV-B had a length of 4629 nt (excluding the 3' terminal polyA tail), which is one nt shorter than that of the AILV-V reported in GenBank. A comparison of the RNA2 coding region sequences of all the isolates showed that AILV-V was the most divergent isolate, with the lowest sequence identities of 83.2% at the nucleotide level and 84.7% at the amino acid level. Putative intra-species sequence recombination sites were predicted among the AILV isolates, mainly involving the genomes of AILV-V, AILV-C, and AILV-B. This study adds insights into the variability of AILV and the occurrence of recombination that may condition plant infection.

SARS-CoV-2 RNA-Dependent RNA Polymerase Follows Asynchronous Translocation Pathway for Viral Transcription and Replication.

Translocation is one essential step for the SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) to exert viral replication and transcription. Although cryo-EM structures of SARS-CoV-2 RdRp are available, the molecular mechanisms of dynamic translocation remain elusive. Herein, we constructed a Markov state model based on extensive molecular dynamics simulations to elucidate the translocation dynamics of the SARS-CoV-2 RdRp. We identified two intermediates that pinpoint the rate-limiting step of translocation and characterize the asynchronous movement of the template-primer duplex. The 3'-terminal nucleotide in the primer strand lags behind due to the uneven distribution of protein-RNA interactions, while the translocation of the template strand is delayed by the hurdle residue K500. Even so, the two strands share the same "ratchet" to stabilize the polymerase in the post-translocation state, suggesting a Brownian-ratchet model. Overall, our study provides intriguing insights into SARS-CoV-2 replication and transcription, which would open a new avenue for drug discoveries.

Complete genome sequence of a new mitovirus associated with walking iris (Trimezia northiana).

The complete genome sequence of a new member of the family Mitoviridae was obtained from walking iris (Trimezia northiana (Schneev.) Ravenna by high-throughput sequencing. This is the first putative mitovirus identified in a monocotyledonous plant. The new mitovirus was tentatively named "walking iris virus 1" (WIV1). The complete genome of WIV1 is 2,858 nt in length with a single ORF encoding a viral replicase (RdRp). The highest level of amino acid sequence identity was 45% to Beta vulgaris mitovirus 1. In the viral replicase, a conserved protein domain for mitovirus RNA-dependent RNA polymerase and six highly conserved motifs were detected, consistent with other members of the family Mitoviridae. Phylogenetic inferences placed WIV1 among members of the genus Duamitovirus (family Mitoviridae) in a monophyletic clade with other plant mitoviruses. Sequence comparison and phylogenetic analysis support the classification of WIV1 as a new member of the genus Duamitovirus (family Mitoviridae).

High-resolution structure of stem-loop 4 from the 5'-UTR of SARS-CoV-2 solved by solution state NMR.

We present the high-resolution structure of stem-loop 4 of the 5'-untranslated region (5_SL4) of the severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) genome solved by solution state nuclear magnetic resonance spectroscopy. 5_SL4 adopts an extended rod-like structure with a single flexible looped-out nucleotide and two mixed tandem mismatches, each composed of a G•U wobble base pair and a pyrimidine•pyrimidine mismatch, which are incorporated into the stem-loop structure. Both the tandem mismatches and the looped-out residue destabilize the stem-loop structure locally. Their distribution along the 5_SL4 stem-loop suggests a role of these non-canonical elements in retaining functionally important structural plasticity in particular with regard to the accessibility of the start codon of an upstream open reading frame located in the RNA's apical loop. The apical loop-although mostly flexible-harbors residual structural features suggesting an additional role in molecular recognition processes. 5_SL4 is highly conserved among the different variants of SARS-CoV-2 and can be targeted by small molecule ligands, which it binds with intermediate affinity in the vicinity of the non-canonical elements within the stem-loop structure.

Atomistic structure of the SARS-CoV-2 pseudoknot in solution from SAXS-driven molecular dynamics.

SARS-CoV-2 depends on -1 programmed ribosomal frameshifting (-1 PRF) to express proteins essential for its replication. The RNA pseudoknot stimulating -1 PRF is thus an attractive drug target. However, the structural models of this pseudoknot obtained from cryo-EM and crystallography differ in some important features, leaving the pseudoknot structure unclear. We measured the solution structure of the pseudoknot using small-angle X-ray scattering (SAXS). The measured profile did not agree with profiles computed from the previously solved structures. Beginning with each of these solved structures, we used the SAXS data to direct all atom molecular dynamics (MD) simulations to improve the agreement in profiles. In all cases, this refinement resulted in a bent conformation that more closely resembled the cryo-EM structures than the crystal structure. Applying the same approach to a point mutant abolishing -1 PRF revealed a notably more bent structure with reoriented helices. This work clarifies the dynamic structures of the SARS-CoV-2 pseudoknot in solution.